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t Spatial wavefront evaluation by intensity relationship. 



A method, and apparatus for accomplishing the method, for determining a phase difference of a 
wavefron: at a first ipupii) plane (P1). the wavefront propagating from the firs: plane to a second i image) 
Diane (P2). The methoc includes the steps of providing an intensity of the wavefront at the first plane . 
measuring an ntensity of :he wavefront at the second plane . and determining the phase difference of 
:he wavefront at the first plane :n accorcance with a transfer T uncttcn that employs the provided 
ntensity of :he wavefront at ;he first plane and the measured intensity of the wavefront at the second 
Diane. A single expression descr.bing an Optical Transfer Function is developed and is shown to involve 
only the unknown aperture phase and known quantities. A solution to tnis expression, achieved hy a 
polynomial expansion technique or by a sampling tecnnique. is shown :o vieid tne ohase at the aperture 
anc. together with the intensity a: the aperture, to defne the aperture wavefront. 
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FIELD OF THE INVENTION: 



T his insertion relates generally to the evaluation of a wavefront of a propagamg sigrai. 
BACKGROUND OF THE INVENTION. 

I: is frequently desirab.e to determine the characteristics of a spatiai pnase profile o*' a wavefron;. Ir optics, 
th s phase (divided by 2 ; is referred :o as an Optical Patr D fference *0 3 D) rrac or aberration -"uncticn. A 
[• rowlecge of the O p D is essential to determine how the wavefront wnl propagate or how well the wave f ront 
can Leimagec. However, in that phase cannot be measured directly, a problem of phase retrieval ;s presented. 

A conventional phase retrieval method employs interferometry, in which two coherent wavefrcnts are 
brought together ( nterered) and the resulting high frequency intensity variatiors. referred to as fringes, are 
-analyzed :c determine t he pnase. A similar interference effect occurs in propagation and imaging. Interference 
effects in images also produce fringes which are referred to as diffraction Tinges or sidelobes. 

Several of t he better K nown properties of mterferometry include t he following. First, only pnase differences 
can be determined and not the absolute phase. Second, for a single static mterf erogram. a sign ambiguity e <- 
tsts. 

Some esser known, but equahy important properties of interferometry include the following. First, omy 
intensities are measured. Phase cannot be measured direct.y, wnich s a result of a funcamental tenet of quan- 
tum mechanics regarding "observables". The phase is computed from intensity patterns or fringes which are 
recorded on fiim or oy optica detectors Second, anc as a result of the immediately preceding property, it s 
an implicit assumption of nterferometr y that the amplitude of the two beams is constant, or at leas; slow y 
varying compared to the fringes, or is explicitly ^ncwr sc that the phase can be computed. This is required so 
that intensity variations, due ;o phase, can be d ( fere r dialed from intensity var a bans due to ampitude. Similar 
i. ufistramts occur in micro wa\, e and FM systems. Third, adequate samel ing of t fie fringes requires at least two 
rumples per cycle 

An overview is now presented of the basic types of interferometers that provide phase retr eva! from im- 
ages. 

An early type of interferometer is tne Micraeison. in which an input beam is split into two parts, one of 
which is a reference and tne other of which samples a surface under test. The two oeams are recombined to 
form an interferogram which allows the phase difference between the surface unoer test and the reference 
beam to be computed. Multiple beam interferometers may be usee .nstead of a Michaelscn-type to obtain 
sharper fringes. However, n all cases tne general approach .s the same: the beam (or beams) from tne surface 
under test is compared to a reference beam (or beams) and the phase is inferred from the resulting fringes. 

A problem that arises during the use of a Michael son (or similar interferometer" is a requirement for high 
stabi ity between the two optical paths and a requirement that no unknown errors exist in the interferometer 
optics or reference. 

Another common type of interferometer is the shearing interferometer, wherein the reference is the sur- 
face under test which is snfted in position The sh ft mav be generated by a shear elate or by other means A 
sign if leant difference h-^twee-i the ^he^p^n n'prfprr.|ijotop mo the Michaels on interferometer is that the sur 
*3c.e is compared wit h -tser i instead of a reference) and spatial phase rj'f fere noes are obtained For sma I shifts. 
\ he phase shifts from a shearing interferometer are summed tr reconstruct the original wavefront phase. Th s 
tvee of nterferometer s less sensitive :c meohaniea stability than ts the Miohaeison. but places mere con- 
straints on the quality and coherence of the source beam. 

Alt rough lot usual y referred tc as nterferometr y anot ner form of interferometry ts imag ng or. more gen- 
erally . propagaticn. Propagation processes are governed by tne laws of diffract on. Although diffraction and 
nterferometrv arc rot generally as see ated. it has been found that c if fraction resu ts f rom the interference cf 
■ cht 'scattered" from the edges of an aperture with the iiqht ransmirteo tnrouqh the aperture. T he resuitinq 
mer fere nee fnnqes. usual v called diffraction siceloces or fringes, are a 'unction of the aperture functicn. sucn 
as sice anc srape. and the distance of the frnqes f rcir :he aperture as we I as a c erratic n .n the aperture. In 
propagating a wa\e*ront. light from a*! carts of 'he aperture s summed together (interfered; to ca cu ate the 
ore page tec wavefront. This :s similar to a shear, ng inter fero meter .n wnich the wave'ront is mer'erec witn user, 
as displaced. Under certa n circumstances the phase in the aperture can be estimated from a measurement 
of an image. 

C orventicna wavefront pnase -ernevai methods from propagation or ■macinc are iterative, successive ap- 
proximation approaches. ~>.ese ap croac-es ar^ generacv n accurate v^c ma, net '-eressar iv converge tc a 
ccrrect solution cue to the e>istence cf sec.or.aarv maxima Vvhen these conventional metre c s do oc^ verge. 
C"^n:p!ice s 'reguenpv . erv n nter*v r ometr. . <■•>■]" S cir 1 : amounts o* "arc ware mav re required \c 
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acccmp isr wavefrcnt cnase retrieval. The provision of acciticnai hare ware is especially disadvantageous m 
airoorne and spaceborne applications, wnere size and weight are important factors. 

I: is thus an ocjec c f the invention to empioy a measLremeni of an image to determine a oer rations, the 5 
ei rnmat ng a reed for :-e hardware associated witn nter f eronetnc apparatus and the staoiiity requirements 
cf maintaining optica; alignments to a Taction c* a wavelength 

I: is a furt her ob ect of the invertion to provide met hed and apparatus for accomplishing Spatial WaveTont 
Evaluation by Intensity Relationships (SWE3IR;. 

SUMMA RY OF THE f NVENTjON 

The foregoing and other problems are overcome and the objects of the irveotion are realized oy method 
and apparatus for determining a phase of a wavefront from a measurement of intensity. The method operates 
in a rapid and exact manner without requiring the use of approximations. The teaching of the invention mav 
be apphec to a large number of applications inducing, but not limited tc, optical, radar, sonar, e-ectron micro- 
scopy, and x-ray systems The method of the invention is rigcrous. exact, fast, anc any sign ambiguities are 
easily reso:ved. Given correct numerical inputs, the method yie'ds results accurate to the computational ac- 
curacy of a data processor employed to perform the methed. 

More specifically the invention teaches a method, anc apparatus for accomplishing the method, for de- 
termining s phase profne of a wavefrcnt at a first plane, the wavefront propagating from the first plane to a 
second Plane.. The met hoc includes the steps of provid in g an intensity of the wavefront at the first plane: meas- 
uring an Intensity cf the wavefront at the second plane; and determining the pnase difference of the waveTont 
at the first plane in accordarce with a transfer function that emp.oys the provided intensity of the waveTont 
at the first plane and the measured intensity of the wavefront at the second plane. 

A snyle expression describing an Optical Transfer Function is developed and is shown to invoke only (a) 
unknown aperture phase anc (b) known quantities. The solution to this expression, achieved by a polynomial 
expansion technique or oy a sampling technique, is shewn to yield the phase at the aperture and, together 
with the ntensity at the aperture, to define the aperture wavefront. 



BRIEF D ESC PIP T I ON OF THE DRAW ING 

The above set forth and other features of the invention are made more apparent in the ensuing Detailed 
Description of the Invention when read in conjunction with the attached Drawing, wherein: 

Figs 1 a- 1 c illustrate an e <ample of diffraction interferometry wherein Fig. "! a shows a relationship between 
a pupil plane and an image plane. Fig. " b shows the pupr phase and intensity , and Fig. 1c snows a dif 
fraction pattern formed at the image plane w;th and without the insertion of a phase pi ate at the pupil plane. 
Fig. 2a shows a relationship due to prooagaticn between a pupil wavefront and an image wavefrcnt; 
Fig. 2b shows a relationship between wavefronts, intensities, transfer unctions, and other optical quanti- 
ties: 

Fig. > i lustrates the relation of the OTF (Optical Transfer Function) to a shearing interferogram: 

F'O. d is a flowchart illu strati no the stecs of a no iv normal pxnsn^inn ;prhr K ^ p ,-f t he ,n vent ion f cr solving 

f or a phase profile o* a pupil from known intensities of the image and pupil planes: 

Fig. h i lustrates an application of a one-dimensional method for solving for a phase profile of a pupil to a 
two-dimensional aperture; 

Fig. na illustrates the operation of a digital sampling approach to solving for the pupi phase prof le 'or a 
five element arrav; 

Fig. nb shows an e>act solution for the five element array 'or ocd and even coefficients 
Fig. " diagrammatically dep cts image samclmg requirements; 

F*q. is a Dlock diagram Hustratinq a closed loop defcrrnable mirror contra s , stem that empi ovs the SWE- 
3IR technique cf the invention; and 

F'Q. 1 ; shows an example of a phase map determined from an intensity map f or a rounc aperture emoedded 
in a square grid. 

DE TA I LE D DE 5 0 RIPT'ON O FT H E INVENT I ON 

O ac :crdance w it h • he riven t:cn. arc o\ wa\ of ntrcducPG". tca rir eadA p e s r c w n t a t c . f 'ra cpc n s ac- 
tually a *Orm of in terercmetr * . Reference s mace tc F'g. 1a. .vnicn :. lustrates an experiment wnereir :ne n- 
tenscv n an aperture \ t pup'i s kept constant cut t he p r ":se s charge c ca e r p a r t o* the a pe r :u re " P\ the 

of a naif-wave : '\v^ piate 2 • f c Ipc As ca r " : r <fe- - F O. verv d if fer-mt c: : f fra d:c n patterns result 
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for the case w incut the pnase pi ate 2 oA; ana witn trie pnase pi ate 2 <B;. 

T ne invention s Dased upon a realization t "at. s;~ce erferem pnase profiles (having the same intensity) 
result in different diffraction patterns, a reversal of this phenomenon enables the phase tc be inferred Tom 
tne diffract on paler n. T ne technique of the invention is thus aptly referred to as Spatial Wave f ront Evaluation 
Dy Intensity Relationships 'S7VE5IR;. 

It is no tec mat SWEBIR reoeires no additional nard ware tc o reduce interference, in that the interference try- 
results from tr-e propagation itself. Furthermore. 3WEBIF is the less sersitive than the conventional nte r - 
ferometrio teenntaues reierrec to apo^e in so far as mechanical stability is concerned, in that SWE3IF: does 
not reouire acditiarai optical components and thus coes not introduce addihona! possiole error sources. Fur- 
thermore, source ccnerer.ee recuirements are minimized and are equivalent to tnat of an equal path interfer- 
ometer. Source quality is similar to that required ir a convent onai shearing interferometer. 

Having thus ntroduced tne SWEBIR technique, a detailec description of same is now provided. 

As was previously discussed, propagation or imaging can oe regarded as a type of interferometry. In con- 
ventional interferometers, a relatively simple algorithm ;s --available for computing phase Torn intensity, ^hat 
is. assuming constant intensity, the intensity (I) is given D\ : 

I = I, [1 + cos <0 - Qv>]. 

From this relaticnsnc the cosine of the phase ditterence can be computed and, therefore, the phase dif- 
ference (with a l- amb gu.ty ) in a non-cenvent onal propagation or imaging type of dif traction irterferometr /, 
the relationship between intensity and phase is mere complicated. However, a def ined relationship does exist. 

The SWE3IR technique of tne invention employs two pas:c assumptions: 

Lthepupiland image intensities are related to one a not ner In a manner tnat involves phase, as was shown 
in F g . 1 ; and 

2. this relationship can be inverted to determine the chase aberratiors in the pupil from a knowledge of 
the pupil anc image intensities. 

To more fully describe the SWEBIR teel muue there is r irst defined the relatiunsh iu between pupil phase 
and image intensity. Ther. there are describee presently preferred methods fur employing these intensities 
sc as to determine the phase and, hence, the aberrations of the pupil. 

Pupil-Image F elatio.nshios 

In Fig. 2a there is snowr a propagation from a first plane (P1 ) to a second plane (P2). It is assumed for 
this discussion that Pi s a pupil and that ^2 is an image, but in general P1 and P2 can represent any two 
planes. The mathematical relationships are shown in Fig. 2b. The two wavefrcnts, at Pi and at P2, are related 
by a Fourier Transform (F. T. ). Mere generally, the two wavefronts car also be shown to oe related by a Fresnel 
Transform if F1. D 2 are not a pupil/image pair. These transforms describe well-known mathematical relation- 
ships that have a unique inverse. Thus, if the wavefront (amplitude and phase) is known in any plane, the wa- 
vefront in the other piane can be uniquely determined. Bv example, the F.T. converts phase information in one 
plane into intensity anc phase variations in the other plane. 

In converting from wavefronts to intensities :n a given plane phase information in that Diane is lost, but 
not the phase nriuced intensity c if fe re noes from t ot he r p: ane Tip ntens tv in t he pupil tstheApertu re Func- 
tion (AFc which is related to *he size ane geometric snape of the pupil The intensity oattern ir the image s 
referred to as a Point Snreae Function (PS Ft (i.e. an image of a o o i o t source spread by diffraction and aber- 
rations in tne pupil). The :ntegra. of the intensity gives the total power Po. which must be equal n both P1 
and P2 so as not to vioate the law of the conservation of energy 

I* the inverse Fourier transform of t he T3F (image intensity ) is taken, the result is referred to as an Optical 
Transfer Function ( OTFi. ~hts function can be comoie>. The magnitude of the OTF is referred to as the Mod- 
ulation Transfer Function (MTF'h The M X F is commork employee in imaging svstems to specK a spatia fre- 
cuencv spectrum passed nv the optical svstern. It is to be noted that bcth the AF and the OTF are defined .n 
the pup I plane ».P1 ) "oq ether, the AF and tne OTF cive two eouations for saving for two unknowns, that is, 
the cup i amphtuae mo the pup 1 pnase. ^hese eouations are. 

Pupil Intensity = I - A.-, and '1) 
O t F = F.'tK - A - e - " * A.e ; , (2: 
where I- and I are the eupn ane mage intensites. respective' v 

It should Le "voted tnat to soive these two eocations eeth h and I must he known, ^his may appear tc 
chfer Tom convention a 1 i r terser ametrv . v. here anu one ntensitv >s measured, however, as was pointed out 
aeove, ;nterferometr > mpiic:tiy assumes t^at f-e lum:ratirc 'ntensdv :s u niform. ~hus. * t s assumed that 
!• s known .-\ ener: tc he a constant over the arerture -,j n i \, tne m a- :~ 'Ctensity I need be measured, ^he as- 
sumption o' constant ; icm mVpc" s reascname *'c m- . -: o :r -t- .r . en . ironments and 'or space i -'sir o men ts 
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oDserv-nc stars. r or the more genera; case, *vhere t- is not constant. oath I- anc b are measurec. 

T he soiuhcn to Equation *1 : ;s relatively sample. A- is taken to be the positive square root of I Eouat.on 
i2; is seen to oe an autocorrelation function. That :S. the aperture wavefrcnt is shiftec. ccmpiex conjugated, 
ccmcarec to itself, sic .ntegratec:. As snown ; n F;j. 3. tms can oe considereo to oe a type of shearing inter- 
fere metr,. ^ne C^F s essential y toe integrate J cut pot of the s near ng interferometer as a function of s near. 
If Equat on ' * > is substituted mtc Equation < 2 ; a smgie ecu at ion resu ts: 



OTF 



= F.T.(I ) = 4T^el + l © ^iT 1 e or 



OTF(r) ^S'Vl^^l^r'- t, e ^*. (r#) -* (r '- r)3 . (3; 



As can he seen. Equation (3) involves only ia; the unknown aperture phase and ih) known quantities. 

The solution to Equation (3) gives the aperture phase 01 and. together with \ i . cefnes the pup i wave- 
front. SWE3IR thus involves methods for solving Equation (3) for the phase in he oupil. based upon the as- 
sumed and/or measured values of h and I-. 



Computational approaches: 



In acccrdance with an aspect of the invention tnere are provided three exact, non- terative approaches to 
solving Equation (3). A *'irst approach is an analytic approach, wherein the phase is represented as an analytc 
function, it being realized that a pO:ynorniai function can be related to common optical aberrations if Zernike 
polynomials are employed. A second approach is a digital approach wherein tne phase is sampled at N points 
and the OTF is sampled at the same points. A third approach is suggested which genera izes tne two above 
mentioned approaches in a systematic: manner. 



Polynomial apprcacn: 



Fig. 4 .s f ow diagram illustrating the polynomial approach for a case where N=6. where the aperture ,s 
one- dimension at (slit), anc where I . = 1 . Starting with t he given expression *or phase in 3 lock A. 0 is expanded 
as a polynomial expansion of order N-1 with coefficients a,. There are N coefficients. a fl ...a*, -. which give N- 
1 powers, plus a reference piston, or cc nstant. phase a,-,. In taking phase differences i Block B) a 0 cancels, leav- 
ing only N-1 coe*f cients. Ths expansion is inserted into Equation .;3). Next, both sices of Equat.on (3: are 
also expanded as a power series for a polynomial expansion) of orcer N-1 and the ccefficients are equated 
{Blocks C and D) T his proviaes N-1 ecuations for the N-1 wavefront coefficients a.. These ecuatiors are then 
manipulated and solved (Block E). After manipulation there are optamed two groups o f equations for the odd 
and even coeff ic ents. ^he odd coefficients are represented bv linear equations ana the even coefficients ty 
q users tic ecuaticn^, as iO iows. 



ODD Cll a. + c 13 a 3 + c 15 a 5 + . . . (4a) 
C 3i a i + C 33 3 3 + C 35 a 5 + 



EVEN c 222 a^ - c 224 a 2 a 4 - c^a^ + ... (4b) 

C 422 a 2 + C 424 a 2 a 4 * C 444 3 4 + ' " ' 

A s seen :v?iow tor on example solution sot for N = h. the linear equations are solved for the ouc ccc 
?nts ana '. quadratic equations ore soiveo for the even coefficients. 
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Odd: 



1/2 ImD 
-1/2 ImD 2 

-1/4 IraD 4 -l/2ReD 3 ImD 2 +l/4 (IroD 3 ) 




a 



a 



1 



4- 3 



+ a 





"2 



1 



a i = E 2' E 3 



a 2 = 2E 3 -5/2(E 2 -E 1 ) 



a 



3 



3/2(E 2 -E 1 ) -E 3 



Even: 



• 4f 1 3ReD . T^h-D-H = F. = + 4a : a, + 4a- 
1 -[ • 1.2ReD : - 1/4(lmD.)-- 4-'5a r - 1ba,a^7 - 1fca-,/9] = = a-. 3 + 4a a.. 5 + 4a-,/7 
a 2 = i 3 4[ - F. i 7/3 • ^ JF. - F- >] a. - 1/2[^ F - - a 2 J (a : + 2a 4 r - F- 

l: can readily be seen from Equation (4b) that since al! terms are quadratic forms, there is a sign ambiguity 
for the f rst even order term. "~his is the same result as in interferometry wherein .t cannot be determined what 
siae of focus tie mterferogram represents, because it is not known if the fringes are increasing or decreasing. 
Tc establish the sign of the fringes, an arrow is frequently placed on the mterferogram or a sign convention s 
adopted. However, once the sign of the first ever coefficient is determined, the s gns of all the other even 
coefficients are unique y determined by the cross terms. 

I: is noted that an expansion in a polysomia of degree N- 1 actually has N soeff cients. a 0 --a n However, 
because of the phase differences. a 0 cancels out. Therefore, tne reference piston pnase nay be chosen arbi- 
trarily (usually se: equal to zero). The one-dimensional approach illustrated in the flowchart of F g. 4 may also 
be applied to twc-dimensions to determine the wsver'ront on an N- 1 by N- 1 grid, as depictec in Fig. 5. To ac- 
complish the method ir: two-dimensions the aperture is masked to pass only one row at a time. Each row has 
N- 1 data points and is represented oy a <N- 1}-" degree polynomial. Tne result of applying the method is to de- 
termine N-1 coefficients for eacn row. There is alsu an arbitrary piston shase coeff cient for each row. that s 
undetermined. N-2 of the N-1 piston coefficients are determined by applying the method to a cclumn. The 
single undetermined coefficient is the piston phase *or the system. Thus, the one dimensional solution illu- 
strated above may be applied to a gnc of 25 points tc determine 24 ohase differences. 

A more general two-dimensional solution does not require masking of the pupil, and is discussed below 
in the subsection ertitied "general approach". 

In practice, tne expansion of the OTF may be accomplished by several techniques. For any given roise 
spectrum there is an optimum polynomial expansion o + ' Equation 3. and it is the derivatives cf these poiynomia s 
that are used by the polynomial approach Other techniques are possible which do net explicitly invove der- 
vatives. but the resulting equation sets (Eauatons <-a) and i4b)i are similar. A noise analysis ndicates that 
tne noise sensitivity in SWEBIR and tne noise sensitvitv ir irrerferometr v are essentia:!--, ecui-alent. 

Sampling approach: 

The sampling approach operates to sample t he wavefront at N points. At earn p nin: the wavefront sample 
is Ae- This sample is inserted into Equation i3i as depicted n Fig. 6a for a *ive element arrav. wherein a = 
uniform amplitude). The results are depicted in Fig nb. wnirh illustrates equat ons o* t he for m (S = shrt. A 



ODD < ASYMMETRIC) terms, set 



0- - 0 _ - V.. 



EVEN , SYMMETRIC} terms, se; 



'0. + 0 . ' 2 - -C. ■ + 0 . . ■ i 2 = y. <5b 
ne ::rHse differences invoive oniy antisvmmetro -oca : acerrations. white : ^e phase averages nvo ve on \ 
^er\; e;e^ aper wooes . It srouti re :. na: ;~e ^ve r - terras :,n;> appear as cos.nes wn cr, nave a 
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sign amoiguity. wnile the ode terms appear only ir exponential form. These VI equations are solved to give 
VI prase values, pus one undetermined reference prase (piston Clear 'y. ' pan pe seen * pat the sampling 
approach exh bits the same properties as tne poivnomial approach. An example solution for f\-2 :s -given ir 
F;g. 6b. If applied to two-d mens ons tr is yie ds 24 pnases. pi us one reference, on a 5 oy 5 grid -25 points). In 
this ac preach tne uniform samp ing f Aes the sensitivity :c noise. To opt;m ze *cr an arDitrary nc se spec:' um. 
nonuniform sampling is employed. 

General apprcacn: 

Each approach discussed above can be e.dended to an arDitrary number o*' points ir one or two dimen- 
sions, however, and as was previous / mentioned, the general two-dimensional approach does not require 
that pupil masking be used tor twe-dimensians. The general approach solves N equations in V unhnowns. 
rather that the simpler case of \ sets of N equators in N unknowns. In all cases the OTF is a -'unction of 
aperture shape. 

It is desirable to pre vide as genera a technique as possible. However, wh lie the technique may be genera: , 
the specif ics cepend upon the a: erture shape noise spectrum, and other constraints (i.e.. computational con- 
straints ). A common example of a computational constraint is the use of the F^T /Fas: Fourier Transform) al- 
gorithm, wmch fixes the grid and mtrocuces aliasing As another example, the geometry of a selected detector 
grid also fixes sampling. 

Given an aperture shape, noise spectrum and sampling recu:rements, an optimal set of orthogonal func- 
tions is determined in which to expand the 0~F. Other sets of functions are then determined in which to expand 
the phase or wavefront. Tnese 'unctions are desired to be near optimal, computational^ efficient, and easi y 
extendable to large N. Because of the non-linear aspect of the problem, a general approach tnat is valid for 
al N is selected. The approach \r. to determine a set of functions which have computational y efficient prop- 
erties. By example, a phase expansion that is orthugona over the aperture, and which retains its ortnogo lal y 
in the OTF equation (Equation w<i). is desirable for computational efficiency. It should be nuteu tnat i: is also 
desirable to fcrmulate a generated set of equations having optimal noise immunity. 

Sampling requirements:: 

The SWE3IR technique requires that the image intensity be measured. In accordance with the Nyguist 
Theorem, at leas: two samples are required per cycle. Two samples per c>c!e in the pupil, wnen transformed 
into image space, requires two samples per - /D. as depicted ir Fig. 7. The widt h of the central image icbe (A) 
is approximately >jD (3dB widt hi. However, many image systems undersampie the image to save detectors. 
In some instances it is still pcssiole to compute the wavefront from an undersarrpied image, but the problem 
is always tnat undersampnnc causes aliasing. ,e., any residual higher order aberrations, such as coma and 
spherical, may contribute :o the computed lower order aberrations (tilt anc focus), 'herefore. auecuate sam- 
pl ng is desiraple, 

Reference s now made to F-g. 8 which llustrates a controller tC for a deformable mirror 12 Defcrmab^e 
m rrrr 12 includes a base 14 and a plurality r.i mirmr sprtmn^ ^6 M'rrrr sections 1 6 are con rented to the rase 
1- hy a plurality of actuators 1S. such as voice cons A pnmt source 20 provides an output *o a beamsplitter 22 
wmch directs radiation Tom the point source 20 to the deformable mirror 1 2. The deformable mirror 1 2 net ines 
the pupil, or aperture, for the system. Raciat on reflecting hack from the de T 'ormacle mirror 1 2 passes through 
the beamsplitter 22 to a magnifying optic 24 and an optional spectra filter 26. A radiation sensitve surface of 
a selector 28 defines the image mare of the system. By example, the detector 28 includes a CCD camera. 
An output of tne detector 28 is provided to an actuator control system 30. the actuator control system 30 in- 
cluding a data processor 32 for executing one of the SWEBIR methocs of tne invention. For the illustrated 
embodiment the SWEBIR method is Me sari pi nq method descrioed m detail aoove. An output of tne control 
svstem J O 'S provided to actuator electronics 34 ,vmcn nas an output 'or driving pe oeformaoie mirror actuators 
\8 As can be seen a closed icoo control svstem is provided wherein the SWEBIR met hoc determines the 
phase, anc hence the a per rations, at the pupil p'ane. In accordance with the ceter mined aberrations * he con- 
tm ^vstem 30 may drive the actuator electronics 34 so a? to null out, or nnimue. the pupii a Per rations. It s 
also witnm the scope of the invention to drve the actuators 18 so as to nduce a desired a per ration at the 
pupil plane 

In operation, a V 1 mage of *>e ro ! rt source 30 ;S magnified -ind measured o\. the detector 2S to accuire 
a measurement of t r e D SF. "^r e maqn r , cation -s crc se n sucn that t he maoe :S sampled at least at tne '-Jvduist 
rate, that is. at a spatiaf ; reoue n cv c' 2 V . A/nere \ ;s the s\ stem '-num-cer I n accordance wit * t he -nventior , 
t detector 2" iTieasures th- rtersr.'. •: ' n e ram 1 piare. Next, one o* tne ■ntens;U patterns is trc-nsformed 
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n:o t h e plane cf the othe r . 5 , exarrpe. :ne point ^mage :s Fourier transformed to y ; e i d the 0~F. wnicr has 
oee" snown to ce the autc-corre.atic n cf tne complex pur i amp ituce witn itseif. Measurement or inference 
of puni intensity prcvijes tne amci tuce. Equation i 3) is tnen solved, as descnoec aDove, to provide the pnase 
oistnout on at tne pupil 

For tms empoaimenr, tne image at tne oupil plane is Kno.vr a cnon. If tne ntersity at tne pupil piane s 
not known t <s witnm the scope of :he invention to provide means, sucn as a Peamspi tter 36 and a separate 
detector 38. for measuring t he wavefrcnt mte is it y at the plane of the pup I . The output cf detector 38 is provided 
to the control system 33 'or use by the SWEBIF'. method, as described aoove. 

Referring to Fig. 9, using a " X 7 square pixel grid as an output of tne detector 28 anc a circular aperture 
having a rad.us of 7, the intensity arounc tne edges of the aperture is tapered accorcmc to a fraction c/ the 
CCD oi.* ei contained wth.n the aperture Of the 49 total pixel?;. 45 of the pixels detect some energy and may 
thus :>e employed to determine tne aperture ohases. This arrangement ideally enables 4 5 Zernike poly norma s 
:o be calculated, It is notea that this is more than the number commonly fitted to interferometr c cata. Fig. ? 
;a-so cepicts an intensity map corresponding to the 45 pixels that receive energy anc further depicts the phase 
mao of the pupil that ts come u ted Tom tne n tensity map. In a computer simulation, tne computed SWEBIR 
phase map was found to equa- the input phase map to seven decimal places. When tne input phases were 
n creased by a factor o: -".the same agreement between the computed and -nput pupil phase maps was ach- 
■e/ec moduio Zpi. The computational time was found to pe ; ess than 0.1 second. 

It should be realized that various mod if ications to the forgoing teaching may be accompl ished. By example, 
it s within tne scope o* tne invention to employ phase masks in the re-imaged pupil to improve accuracy for 
near-perfect optcs. The phase masks change the "interferocram" Tom a null fringe to a til tec or displaced 
fringe, providing optimum- sensitvity. Various such masks are possible. Examples include, but are not limited 
tp. tilts, aefocus. and corner pi*els shifted bv 90". 

It is also within the scope o*' the invention to employ variuus noise reduction techniques. These include, 
nut are not limited to, a tec hnique of ado i no small errors to tne computational met hods, linearizing tne equations 
; i the errors, and solving the resulting linear equations by regression techniques which minimize the uu:st*. 
^his is preferaolv accomrris^ed after each row or column computaron. sucn as that illustrated in Fig. 5, bet 
■'■an also be accomplished at the end of the computation To prevent normallv computaole numbers from be- 
ccming ncomoutable due to the presence o" noise (i.e. cosine > unity), ar error bar is employed sucn that: 
■'' zosme < 1 use cosine 

*+k cosine > 1 cosine = 1 

cosine > 1 +k cosine out of range. 

In summary, it has been shown that the phase of a wavefront at the pup i can be reconstructed from the 
■mage, provided the puoii intern: it y prof ile is known. As an example, a one -dimensional pro Diem with a uniform 
aperture it) is solved by the use cf two teenmques an analvt c polynomial expansion technque and a digital 
sampled technique. Both teenmques have been snown to exhibt the same properties, wnicn are identical to 
those of interferometry . Specif cally. the SWEBIR technique of the invention has been, shown to exhibit the 
'o lowing characteristics in common w;th conventional intererometric measurement techniques. 

i;a) only phase cifferenoes are determined; 

if-) a sign nmniouity ^xi^t^ -or fir^: even order aherrat-on. 

i: ) orvy intensities are measured: 

(dl the intensity profile in t^e aperture is required to be known, and is assumed to be constant or is meas- 
ured: and 

(e) adequate sampling is required, for example, two samp es. cycle. 

It nas also been shown that tne one-cimensionai technique nay be extended to two-dimensions 
Tne methods of the invention nave a nunDer o*' appitcat or s. In optics, by example. SWEBIR can eliminate 
: ne need for interferometers. T hus, one troav determine t he aberrations of a s pace- b a sec telescope oy ooser v - 
nc a star and measuring :ne maqe < p SFc For around-based acpl ications. atmospnenc measurements and 
compensation can t:e acccmpisnea ^ measuring ooth the puoi! anc imaqe ntensites. SWEBIR rrav aiso Pe 
emptcved to measure the ;ur mature o* aspher-cs without a nu lens. With known optics. SWEBIR can be em- 
ployed to construct q scene < n:;udinc onase) wnei Tunicate;] with coherent hqht. One appi cation is tne :n- 
sreovon of m'crocrcu ts + cr defects Cther applications c;c;uoe tne measurement and correction of prase er- 
rors n ngr frequency side iool- inq radars. In general, any measurement using waves, such as optics, radar, 
x-'avs. acoustics, scnar. arc electrons i electron microscopes) in wh on an mtensitv is measured and a chase 
s :o oe determined, mav oe r ~ef iciailv emmm.- the teacmnq of tne '^venticn. 

* * us wm e 'he invention n^s tee" : -r ' cu ar . snown am: described with resoect to a preferred embodi- 
ment thereof, it wri ce m tier- toed bv rcse sk Dec r ar" t^a: c" a-~ >-< n 'enr and deta :s -v.av re made 
tnere'^ wnnou* depart;-^: *ron :^e -cor e r\: sr :r.t ./ t :^e '^v error 
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Claims 

1. A met nod for determ ntng a pras-? prof le of a wavefront at a first pi ane, tne wavefront propagating *rom 
:ne f :rst plane to a seconc Diane, comprising the s'eps of: 

providing an intensity o* the /.avefront at the first piane. 
measuring an ntensty of" the wsvefront at t~e second plane; and 

determining the phase pro* le of the wave Yon t at : he f rst plane ;n accorcanse <vitn a transfer *u no- 
tion that employs t ne provided intensity of the wavefrcnt at the first p ane and the measured i men s ty of 
the wavefront at the second p ane. 

2. A metnod as set *orth m Claim 1 wnerein tne step of providing includes a step of measuring t he intensity 
of the wavefront at the first plane. 

3. A metnod as set fnrtn n C aim 1 and including a step of determinnc the transfer function by taking an 
nverse Fourier transformation of the measured ntensity at the se:ona plane. 

4. A met hot: as set fcrt h in C laim 1 wherein the step o f providing includes a step of determining the intensity 
of the wavefront at an aperture of an apparatus t hat receives the wavefront. and wherein the step of meas- 
uring incudes a step ot measuring tne ntensity of the wave f ront at an image plane of tne apparatus. 

5. A metnod as set forth r Claim 1 wherein the step of determining induces a step of performing a poly- 
nornia expansion, for N-i aberration coefficients, by the steps of: 

ca cuiating U-~ derivatives of the transfer function; 
obtaining <N-1 i/2 linear equations for cdd aberrations, 
obtaining <r-M)/2 Quadratic ecuations for even coefficients; and 
so ving the equabons to determine the phase profie. 

6. A met hoc as set *"crth in Claim 1 wherein the step of determining mc;udes a step of sampling the intensity 
at N points at the second p ane. and further includes the steps of: 

from the sample intensity. 

obtaining <r*-1»/2 equations for odd aberrations, 
obtaining <N-1)/2 equations for even aberrations, 

obtaining ■;N/2''-* equations to resolve sign ambiguities for the even aberrations, and 
so vircj the equabons to determine the phase prof le. 

7. A metnoc as set forth in Claim 1 where: n the step cf measuring includes a step of operating a wavefront 
ntens ty detector to measure the wavefront at a plurality of points 

8. A method as set forth in Claim 1 ^vhere n the step cf measuring nciudes a step of operating a waveYont 
ntens tv detector to measure the wave'Yont at a rate equal to or exceeding the Nvquist rate. 

9. Apparatus far ceterm;r,ng a pnase map of a wavefront at a cupil blare, the wavefront propagating Yom 
the pup I plane to an image plane, comprising. 

first means for measuring an intensity of the waveYont at the image plane, and 
means, having ar input coupled to an output of said measuring means, for determin ng the phase 
difference c*'the wave r Yont at the pupil plane in accordance witn a transfer function thatempioys a known 
or measure': intensity of the wavefront at the pupil plane and the measured intensity of the wavefront at 
the image pane 

10. Apparatus as set forth ! n C aim 9 and fur t her inducing ^cu iJ means fur i red sunny the Ttersitv o' the 
wavefrort at the pupii oiane, said second measurng means "av ng an : u:cu: ccuc ed to sate determining 
means i c r prov ding the measured intensity tnereto. 

11. Apparatus as tort h iM C aim ^ wherein said ietei r : means operates m ac:orianc:e with a deter- 

[Tiinpri optical trr-nsf^r Yj notion obtained m taking an nverse Fourier transfer mat inn nf the measured tt- 
•ensitv at imace mane 

12. Apparatus as se* 'or: w - Ca.m - wherein said determ rmc means ncudes data processor means *c r oe r - 
; ormmo a polynomial expansion. r or V* aperraton coe f f cients. a^d -nc udes 
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means fr.r calculating N-1 derivatives of tne :ransfer function; 
means 'or co taming < N-1 >.-2 linear ecuaticns for oca aoerratiors. 
mears . r o b ta i r i n g < N - 1 ; ■ 2 Quadratic ecuaticns f cr even cc ef f i c ; e n t s ; ana 
means f -::r soivmc t^e ecuaticns to deter nine tne c^ase map. 

Apparatus as set fortn in Claim 9 vvnerem saia cetermining means includes a cata processor means for 
sam.pl ng tne measure intensity at N pcmts at the image plane, and further induces: 

means fr r obtain no i N-1 ) I ecuaticns for odd a Der rations: 

means fc r obtaining < N-1 ).Z equations for even aberrations. 

means f<: r obtaining (N/2) - 1 equations to resolve sgnambiguitesforthe even aberrations: and 
means for solving the equations to determine the phase map. 

Apparatus as set forth in C aim 9 wherein said measuring means inciuaes a wavefront intensity detector 
*'or measuring the wavefront a: a plurality of points. 

Apparatus as set forth in C aim 9 wherein said measuring means includes a wavefront intensity detector 
that measures tne wavefront at a rate equal to or exceed ng tne Nyquist rate. 
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